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Exact  Significance  Testing  with  Biased 
Coin  Randomization 


Myles  Hollander  and  Edsel  Pena 


■assa. 


Abstract 


For  a  clinical  trial  where  two  treatments  have  been  assigned  sequentially 
to  the  patients  via  Efron’s  (1971)  biased  coin  design*  a  recursion  procedure  is 
derived  for  obtaining  the  exact  randomization  distribution  of  a  class  of  test 
statistics.  This  enables  one  to  perform  exact  significance  tests  of  the 
hypothesis  of  no  treatment  difference.  The  randomization  distribution  of  the 
statistic  is  conditional  on  the  imbalance  of  the  treatment  allocation.  It  is 
illustrated  that  if  the  analysis  is  performed  as  if  complete  randomization  was 
used*  conservative  and  anticonservative  errors  can  be  incurred.  The  applica¬ 
bility  of  the  test  to  censored  data  is  also  discussed. 
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1.  INTRODUCTION 


In  clinical  trials  to  compare  the  efficacy  of  treatments  A  and  B,  patients 
typically  arrive  sequentially  and  upon  arrival  are  assigned  a  treatment  or 
otherwise  excluded  from  the  trial.  Several  designs  are  possible  for  assigning 
the  treatments.  Complete  randomization  assigns  treatments  on  the  basis  of  the 
outcomes  of  independent  tosses  of  a  fair  coin,  whereas  systematic  design 
alternately  assigns  the  treatments  after  the  first  patient  has  been  assigned  a 
treatment.  The  complete  randomization  design  is  optimal  in  minimizing  selection 
bias  (see  Blackwell  and  Hodges,  1957),  the  bias  introduced  when  the  experimenter 
can  predict  with  high  probability  the  treatment  an  incoming  patient  will  receive. 
It  also  controls  accidental  bias,  the  effect  of  covariates,  time  trends,  etc.,  on 
the  statistical  analysis  of  the  resulting  data.  However  it  suffers  from  the 
defect  that  it  could  produce  highly  imbalanced  treatment  allocations  especially 
for  small-sample  trials  (Efron  1971  and  Pocock  1979).  Such  an  imbalance  may 
decrease  the  efficiency  of  statistical  procedures,  and  may  also  lead  to  less 
credible  results  (Halpem  and  Brown  1986;  Smith  1984b).  In  contrast,  the 
systematic  design  is  optimal  in  balancing  the  treatment  allocation,  but  is  the 
worst  design  in  controlling  selection  and  accidental  bias. 

Efron  (1971)  introduced  the  biased  coin  design  with  bias  p,  abbreviated 
BCD(p)  for  convenience.  This  design  is  a  compromise  between  the  complete 
randomization  and  systematic  designs.  Let  Tj,  T2,  ...  denote  the  sequence  of 
assignment  variables  with  T.*0  or  1  according  to  whether  the  i**1  patient  receives 
treatment  A  or  B,  respectively.  With  DQ=0,  define  for  i  =  1,  2,  ... 


Under  the  BCD(p)  where  0.5<p<l,  the  probability  that  the  (i+l)*  patient  receiver 
B  is  p,  0.5  or  q=l-p  depending  on  whether  <,  =,  or  >  0,  respectively. 
Wei(1977,  1978)  generalized  this  by  allowing  the  assignment  probabilities  for  the 
(i+l)th  patient  to  depend  both  on  i  and  D.,  Efron  (1971),  Wei  (1977,  1978)  and 
Smith  (1984a,  b)  demonstrated  that  these  designs  compare  quite  well  with  complete 
randomization  in  controlling  selection  and  accidental  bias  in  addition  to  having 
better  balancing  properties. 

We  confine  our  attention  to  a  clinical  trial  where  treatments  A  and  B  have 
been  assigned  to  n  patients  via  the  BCD(p),  and  study  the  randomization  test  of 
Hq,  the  null  hypothesis  of  no  treatment  difference.  We  consider  test  statistics 


of  the  form 


(1.2) 


so  that  the  two-sided  version  of  the  randomization  test  rejects  Hq  whenever 
is  either  "too  small"  or  "too  large."  In  (1.2),  aj,  ...,  an  is  a  nonrandom 
sequence  of  scores  associated  with  the  sequence  of  patient  responses  Xj,  ...,  x^, 
and  these  scores  are  typically  functions  of  the  ranks  of  the  xi's.  In  deciding 
whether  Sn  is  too  small  or  too  large,  Cox(1982)  suggested  taking  the  randomization 
^•^ribution  Sjj  over  those  treatment  allocations  with  the  same  or  nearly  the 
same  terminal  imbalance  as  the  observed  allocation.  In  accordance  with  this 
suggestion,  our  results  are  conditional  on  Dn,  the  terminal  imbalance  of  the 
treatment  allocation. 

In  Section  2  we  utilize  the  Markovian  structure  of  DQ,  Dj,  ...  to  derive  a 

recursion  procedure  for  obtaining  the  exact  distribution  of  S  .  conditional  on 

n 

Dn»m.  This  procedure  enables  one  to  perform  exact  significance  tests  of  HQ. 

In  Section  3  we  present  procedures  for  computing  higher-order  probabilities  of 
®g»  ®l*  •••  suggest  possible  approximations.  These  probabilities  are 


needed  in  the  recursion  of  Section  2.  Section  4  illustrates  the  consequences  of 
ignoring  the  BCD(p)  allocation  and  instead  acting  as  if  complete  randomization  v;a^ 
used;  while  Section  5  considers  the  applicability  of  the  randomization  test  to 
censored  patient  responses. 

Whereas  our  work  focuses  on  exact  results,  previous  research  dealt  mostly 
with  asymptotic  results.  Efron  (1971)  presented  an  asymptotic  argument  showing 
that  if  one  ignores  the  BCD(p)  allocation  and  instead  acts  as  if  complete 
randomization  was  used,  the  randomization  test  of  could  be  conservative  or 
anti conservative.  Halpem  and  Brown  (1986)  concluded  on  the  basis  of  their 
simulation  results  that  in  the  case  where  the  patient  responses  are  binary,  the 
classical  2x2  x  -test  should  not  be  used  to  compute  significance  probabilities 
when  the  observed  responses  exhibit  a  strong  trend.  Smythe  and  Wei  (1983) 
derived  the  asymptotic  null  distribution  of  Sn  under  Wei's  (1977)  urn  design. 

Cox  (1982)  and  Smith  (1984b)  discussed  the  randomization  tests  of  HQ  for  special 
cases  of  Wei's  (1973)  biased  coin  designs.  Wei,  Smythe  and  Smith  (1986)  derived 
the  asymptotic  null  distribution  of  the  k-treatment  version  of  (1.2)  when  the 
treatments  are  assigned  via  a  k-treatment  version  of  Wei’s  biased  coin  design. 

2.  THE  RANDOMIZATION  DISTRIBUTION 

Let  z  denote  the  set  of  integers,  the  set  of  positive  integers,  and  set 
Z^  *  Z^  u  (0).  By  the  defining  property  of  the  BCD(p),  the  process  DQ,  D^,  ..., 
in  (1.1)  is  a  homogeneous  Markov  chain  with  state  space  Z.  It  has  stationary 
transition  probabilities 


PiJ  s  Pr(Di-j|D0»i), 


Ul.  i €l 


if  j  =  *1,  i  =  0 

if  (jai+1,  i<0)  or  (j*i-l,  i>0) 
if  (j*i-l,  i<0)  or  (j-i+1,  i>0) 
otherwise. 


(2.1) 


(P?  ieZ,  j  e  Z}  be  the  n**1  order  transition  probabilities  of  this  chain, 
and  represent  the  randomization  distribution  of  Sn,  given  Dn  =  m,  by 


h^(s) 


|  Pr(Sn=s|D0=0,Dn=m). 


s  €  Pft  m  >  0 
m’  0,m 

otherwise 


where 


n  n 

A  *  (  [  aiti:  €  tO*1*  3X1(1  2  £  t.  -  n  =  m). 

i=l  i=l  1 

Then  we  have  the  following  procedure  for  computing  {hn(s)).  For  continuity 
of  presentation,  proofs  of  the  theorems  are  deferred  to  the  Appendix. 

Theorem  1.  If  m  e  Z  with  pJ  _  >  0  and  s  e  A®  then 

~  u,m  m» 

hm‘s>  * 

Jq 

where  (J^Cs)}  satisfy  the  recursion  equation 

Jmts)  '  Vl  *■» . . 

with  initial  and  boundary  conditions  jj^(s)  *  1  if  m=  0  and  s  =  0 ,  0  otherwise;  and 
where  Ymap,  .5,  q  according  to  whether  m  <,  =,  >0,  respectively. 


An  immediate  consequence  of  this  theorem  is 


I 

I 

I 


\ 


(2.2) 


w^ere  Iu  denotes  summation  over  all  u  e  A^.  Combined  with  Theorem  1  this  shows 
that 


S£<- 

Using  (2.2)  to  obtain  p”  m  requires  knowledge  of  J^(u)  for  all  ueA^.  But 
to  compute  significance  probabilities  one  usually  needs  only  those  values  of 
h”(u)  and  j”(u)  for  u  beyond  the  observed  value  of  Sn.  In  Section  3  we  therefore 
present  a  different  method  for  computing  Pq  which  does  not  require  knowledge 
of  {J>)>. 

The  next  result  allows  us  to  restrict  attention  to  conditional  randomization 

distributions  with  D  =m2  0. 

n 


Theorem  2. 


!f  meZ^  with  Pj^m  >  0  and  s  e  aJJ,  then  hjj(s)  =  h"n( .X  aA-s) . 


Observe  that  when  a^  =  rank(x^)  the  recursion  in  Theorem  1  has  some  similarity 
with  that  of  Mann  and  Whitney  (1947)  for  the  Mann- Whi tney- Wi 1 co xon  U  statistic. 
However,  in  contrast  to  the  distribution  of  U,  (h”(s)}  is  not  invariant  with 
respect  to  permutations  of  the  scores  a^,  ...,  aR.  This  is  illustrated  by  Table  1 
which  summarizes  the  conditional  randomization  distribution  of  S  ,  under  BCD(2/3), 
for  all  permutations  of  the  ranks.  (For  economy  of  space  we  just  list  the 
distributions  for  n*4  and  m=0.) 


Table  1.  Randomization  Distributions  of  S  for  the  24  Rank  Permutations  under 
BCD(2/3)  with  n  *  4  and  m  =  0. 


Rank 

Sequences 

1234 

1243 

2134 

2143 

4321 

3421 

4312 

3412 

1324 

1342 

3124 

3142 

4231 

2431 

4213 

2413 

1423 

1432 

4123 

4132 

3241 

2341 

3214 

2314 

Values  of  S_ 


_2_  _3_  6  3  2 

16  16  TeT  16  16 


We  developed  a  FORTRAN  subroutine  that  implements  the  recursion  procedure 
in  Theorem  1  for  the  case  a^  =  rank(x^)  and  Dn  =  iaa0,  Interested  readers  could 
obtain  this  program  by  writing  to  Gdsel  Pena,  Using  the  Cyber  730  computer,  the 
computer  time  required  by  this  program  to  obtain  the  conditional  randomization 
distribution  of  Sn  for  n  =  6,  12,  18,  24  and  30  with  m=0  and  p=2/3  were  0.15, 
0.81,  3.16,  9.14  and  20.75  CPU  seconds,  respectively.  Note  the  exponential  rate 
of  increase  of  the  time  as  n  increases.  In  a  forthcoming  report,  we  present 
large-sample  approximations  that  enables  one  to  obtain  approximate  p-values  when 
the  recursion  in  Theorem  1  is  not  practically  feasible. 

3.  HIGHER-ORDER  PROBABILITIES 


We  now  present  recursive  methods  for  computing  the  higher-order  transition 
probabilities  of  the  process  DQ,  Dj,  ....  Let  Yj,  Y2,  ...  be  identically 
distributed  and  independent  random  variables  with  Pr(Y,=l)  *q  =  l-PrfY.=-l). 


With  Wq  -  0,  define  *  Yj  ♦  for  ieZ^.  The  process  WQ,  Wj,  ...  is  the 

unsymmetric  random  walk  with  negative  drift.  For  j  <■  Z®  and  n  e  2^,  let 

*0,1  -  Prti=VDi’i):  VJlD0=01' 

b0,J  s  Pr(i2ltWiM>  :  V^V01- 

Furthermore,  let  C(b,a)  denote  the  number  of  combinations  of  a  items  taken  from 
b  items. 

Then  the  higher-order  transition  probabilities  satisfy  the  following: 


Theorem  3.  For  n  e  Z  and  nuZ  . 
'  -  ***♦* 


(i)  p«.o "  j/o.o  • 

(ii)  ^  *  ^K.h  *  j/o.o  pS;! 


Theorem  4. 


(i)  ^0,0  *  0  *0,0  ”  C(2n-l,n-l)pnqn”  for  neZ^, 

(ii)  bo,m  =  U«*l)/(n*l)>  C(n*l, (ntm+2)/2)pCn“m)/2  q(n+m)/2 
for  ne  7?,  me  Z°. 

Efron  (1971)  obtained  the  stationary  distribution  of  |Dq|,  | | . 

From  this  distribution,  the  stationary  distribution  probabilities  of  n^, 
denoted  by  (ir^,  jeZ)  are  found  to  be 

*0  *  Cp-q) /C2p)  and  it  .  »  tt  .  =  (p-q)/(4p2)(q/p);*‘1,  j  e  Z  . 


Since  Dg,  Dj,  ...  has  period  2  it  follows  that  the  limiting  values  of  Pq1^ 
and  pon2m-l  are  an(*  ^ir2m-l*  ^respectively.  For  large  n  we  could  approxircai  ? 
Pq  m  in  Theorem  1  by  these  limiting  values.  Computations  show  that  these 
approximations  are  quite  good  for  sample  sizes  of  at  least  30.  This  is  illus¬ 
trated  by  Table  2  which  summarizes  the  exact  values  of  p!?  for  6sns  30  and 

u,m 

Os  ms  6  under  BCD(2/3).  These  values  were  computed  using  Theorems  3  and  4. 

The  last  row  contains  the  limiting  probabilities.  Notice  the  close  agreement 

30  29 

between  Pn  and  2n  when  m  is  even,  and  P~  and  2 tt  when  n  is  odd. 
u,n  m  ’  0,m  m 

4.  A  SIMULATION  STUDY 

A  computer  simulation  was  performed  to  illustrate  the  consequences  of 
ignoring  the  BCD(p)  allocation  and  instead  acting  as  if  complete  randomization 
was  used.  This  simulation  was  done  on  a  Cyber  730  computer  at  the  Florida  State 
University  Computing  Center.  The  uniform  random  number  generator  used  was  the 
intrinsic  routine  RANF. 

Five  hundred  replicates  were  generated  of  the  following  experiment.  Each 

experiment  consisted  of  generating  n=15  independent  uniform  (0,1)  variates 

Xj,  ...,  x15  and  obtaining  their  associated  ranks  aj,  ...,  a15,  and  then 

generating  the  treatment  assignment  variates  tj,  ...,  tJ5  via  the  BCD(2/3). 

After  stratifying  these  500  replicates  according  to  their  value  of  D1& =  m,  the 

conditional  randomization  distributions  (h^(s)}  were  obtained  using  the  FORTRAN 

m 

program  mentioned  in  Section  2.  For  each  of  these  distributions,  the  a-level 
conservative  critical  values  Sq  corresponding  to  the  one-sided  test  which  rejects 
Hq  when  S^j.  is  small  was  determined.  The  significance  levels  were  set  to 
0.01,  0.05  and  0.10.  Note  that  Sq  is  conservative  in  the  sense  that 
Pr(S15  s  sQ|  DQ  «  0,  D15  =  m)  s  a  and  Pr(Sl5  s  sQ  ♦  lj  DQ  =  0,  D15  =  m)  >  a. 


9 


Table  2.  Exact  Higher-Order  Transition  (P^  )  and  Limiting  (2-rrm)  Probabilities 

of  Dq,  Dj,  ...  under  BCD(2/3). 


m 

n 

0 

1 

2 

3 

4 

5 

6 

6 

0.5597 

0.1893 

0.0288 

0.0021 

7 

0.4060 

0.0823 

0.0110 

8 

0.5413 

0.1902 

0.0347 

0.0041 

9 

0.3975 

0.0866 

0.0143 

10 

0.5300 

0.1902 

0.0384 

0.0058 

11 

0.3913 

0.0890 

0.0167 

12 

0.5224 

0. 189£ 

0.0408 

0.0071 

13 

0.3878 

0.0905 

0.0183 

14 

0.5171 

0.1896 

0.0424 

0.0081 

15 

0.3850 

0.0915 

0.0195 

16 

0.5133 

0.1893 

0.0435 

0.0088 

17 

0.3828 

0.0921 

0.0204 

18 

0.5104 

0.1890 

0.0443 

0.0094 

19 

0.3812 

0.0925 

0.0211 

20 

0.5083 

0.1888 

0.0449 

0.0099 

21 

0.3800 

0.0928 

0.0216 

22 

0.5067 

0.1886 

0.0453 

0.0102 

23 

0.3790 

0.0931 

0.0219 

24 

0.5054 

0.1884 

0.0456 

0.0105 

25 

0.3783 

0.0932 

0.0222 

26 

0.5044 

0.1882 

0.0459 

0.0107 

27 

0.3777 

0.0933 

0.0225 

28 

0.5036 

0.1881 

0.0461 

0.0109 

29 

0.3772 

0.0934 

0.0226 

30 

0.5029 

0.1880 

0.0462 

0.0110 

2it 

m 

0.5000 

0.3750 

0.1875 

0.0937 

0.0468 

0.0234 

0.0117 

Table  3  is  a  summary  of  the  results  of  this  simulation.  The  second  column 


shows  the  number  of  replicates  that  have  *  m.  The  last  three  pairs  of 
columns  show  the  percentages  of  replicates,  for  each  value  of  m,  that  have 
a* level  critical  value  sQ.  We  excluded  the  cases  o  =  -5  and  m  *  S  since  there 
were  only  IS  and  9  replicates  in  each,  respectively.  Those  values  of  Sq  that  are 
superscripted  by  an  asterisk  are  the  a- level  critical  values  under  complete 
randomization. 

Table  3  shows  that  if  one  ignores  the  BCD(p)  allocation  and  instead  performs 
the  randomization  test  employing  the  critical  value  derived  under  complete 
randomization,  the  test  could  be  conservative  or  anti conservative.  Conservatism 
is  illustrated  by  the  cases  m=-l  and  m=l  with  a  =  0.01,  while  the  latter  is 
manifested  by  the  cases  ra=  -1  and  m  =  1  with  a  =  0.05  and  a  =  0.10. 


5.  APPLICABILITY  TO  CENSORED  DATA 

The  simulation  study  in  the  preceding  section  dealt  only  with  uncensored 
responses.  However,  the  test  discussed  here  can  also  accomodate  censored  data, 
albeit  from  a  more  restricted  censorship  model  than  is  typically  assumed.  Let 
Xj,  ...,  Xn  be  the  independent  response  variables  of  the  patients,  and 
V.,  ....,  Y  be  the  sequence  of  independent  censoring  variables,  with  the  tXj) 
independent  of  the  {Y^.  In  the  typical  nonparametric  two-sample  censorship  model 
it  is  assumed  that  Yt  conditional  on  ^  =  0  has  distribution  Gj,  and  Y^  conditional 
on  Ti=  1  has  distribution  G2,  where  Gj  and  G2  are  unspecified.  For  our  randomi¬ 
zation  test  to  be  valid  in  this  censored  situation,  we  need  to  impose  the 

requirement  that  (This  "equal  censoring  distribution”  restriction  may 


Table  3.  Number  of  Replicates  With  =  m,  and  Percentages  of  Replicates  For 
Each  Value  of  m  with  a-Level  Critical  Value  s^. 
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not  be  unreasonable  in  clinical  trials  where  patient  arrival  and  treatment 
assignment  is  sequential.)  Under  this  restriction,  the  test  could  then  make 
use  of  a  wide  variety  of  choices  for  the  a*s.  For  example,  Gehan's  (196S) 


method  of  assigning  the  scores  is  as  follows:  Let  (Z ^,6^),  ...,  (zn»5n)  b® 
the  observed  censored  data  where  =  minCX^Y^)  and  5i =  I(X.^Y.).  For 


i,  j  *  1 . .  define 


0  if  (Z.SZ 6.=1) 


*  |  0.5  if  (Z.SZ..,  «i=0)  or  (Z.>Z y  6^0) 


if  (Z.>Z.,  6.=1). 


The  Gehan  scores  are  obtained  by  letting  a.*l+  \  n. .,  i*  1,  ...,  n. 
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APPENDIX:  PROOFS 


Proof  of  Theorem  1.  Let  me  Z  with  pU  _  >  0  and  s  e  A_.  Then 
■■■  ~  o,m  m 


h"(s)  *  pr(Sn=s,Dn  l=m-l|D0=0,Dn=m)  ♦  Pr (Sn=s , Dn_ ^m+1 1 DQ=0 , Dn«m) 


Pr(Sn-l*s"an»Dn-lBm"1l  V °»Dn=m)  +  Pr(Sn-las»Dn-rm*1l  Dos°»Dn=m) 


*  Pr(Dn-l=ra"1 1 D0*° » Dn=m) Pr(Sn- l*s“a J  D03° » Dn- 1*®"1 *  Dn=m) 


+  Pr(Dn_i'sn,+1  |  D0=0,Dn=m)Pr(Snl=s  |  D0=0,Dnl=m+l,Dn*m) , 


Conditional  on  Dn  S^_1  and  Dn  are  independent  by  the  Markov  property  of 
V  Dl*  ••••  Dn*  pu*hermore»  Pr(Dn-1-ra-l|D0=0,Dn=m)  *  and 

Pr(Dn-l“m+1lD0=0,Dn=Tn)  3  po’Ll(1-Ym*l)/P0,m-  Therefore* 


'\vs\v\%v 


'o.mO*)  *  VlC-lCll<s-V  *  u«*»*  Jm<s>  * 

k  I,  lr 

Pn  h  fs)  we  obtain  the  recursion  equation  for  {J  (s) },  The  initial  and  bounda;/ 
u,m  in  in 

conditions  follow  fTom  the  fact  that  h2(s)  =  1  if  m=0  and  s  =  0,  0  otherwise,  and 

nl 

P~_  =  1  if  m=  0,  0  otherwise, 
u,m 

Before  proving  Theorem  2  we  first  prove  the  following  lemma. 

Lemma  1,  The  process  DQ,  Dj,  ...  is  symmetric  in  the  sense  that 

Pr(i5l(Di=di)l Vd0}  ■ 

for  every  dQ,  ...,  dR  with  d^Z. 

Proof  of  Lemma  1.  From  (2.1)  we  obtain 

Pr(Di.rdi.ilDi-di)  ■  Pr(Vr-diJDi-di)- 

By  the  Markov  property  of  {D^}  we  have 

Pr{i21(Disdi)lDosdo>  *  j1Pr(Disdi,Di-iBdi-i) 

=  i21Pr(Di*’dilDi-i*“di-i)  =  Pr{i21(Di*“dilD0*"do}'  ^ 


Corollary  1.  pJ^  =  P-i,-jf  i  €  n£?£* 


Proof  of  Corollary  1.  Follows  from  Lemma  1. 


Proof  of  Theorem  2.  Given  (t.,...,tn)  let  (t£,...,t')  *  (1-t  .  .ti-t  ).  1,16,1 

n  "  IT  n  r  v 

2  J  t.  -  n  a  m  and  £  a.t.  =  s  if  and  only  if  2  £  tf  -  n  *  -m  and  I  aitf  *  \  a^  -  s. 
i=l  1  i=l  1  i»l  i=l  i«l 


By  Lemma  1  and  Corollary  1  it  follows  that 


Pr^i2l(Ti*ti)  t  Do“°*Dnam^  =  Prti21(V1_ti)*V0'Dn-mK 


n  n  1 

Consequently,  h  (s)  =  h_n(  J  a^-s). 

i=l 


•>  '>  *jVji  »j,  *  •  •  >  •»  '  *  **  *ji  v  >  ■>  •  •  ’j.  ? 
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Proof  of  Theorem  3.  To  prove  (i)  we  have 


n  k-1 


■m  ■  PrCk“i{jiitDj'0);  v°'V0,*v01 

s  I  PrtVfD^O);  D^-OlDQ-OJPrCD^lD^-O) 
k*  1  3 

using  the  Maikov  property  of  {D^}.  Since  is  homogeneous,  (i)  follows. 


meZ.,  we  have 


'O.m  1  Prli£lCDi>0) :  DiHD<f0} 


*  PrSc^l{i2i(Di 5(0)1  Djc*°»Dns,n}lDo“0j 


=  (.5)Pr{  n  (D^l);  D^mlD^l} 

♦  l  Pr{kn|(D.^0);  D^OlD^OjPrCD^mJ^O) 
k=l 

«„.!*»-! I W0-0)  . 

since  Dj,  D2,  ...  is  stochastically  equivalent  to  WQ,  J)^,  ...  when  Di  a  1  for 
i*  1,  2,  ...»  and  DQ,  ...  is  homogeneous.  Recalling  the  definition  of 
b0~m-l  we  obtain  (!*)•  II 

Proof  of  Theorem  4.  That  *0  is  immediate  from  the  fact  that  {D^}  is  of 

period  2.  On  the  otherhand, 

£0?0  *  fCD^l.o^1 . D2„_121'D2n'°IV0) 

*  MDj-1  .Oj*-l , . .  .D2n.^-l,I>2n-0 1 D0.0)  . 

By  Lemma  1  these  probabilities  are  equal,  hence 

fM  *  Pr'D221 . D2n-l21’D2n-°lD0-0’Dl*» 

*  Pr(W.>-l , . . . ,W2n_2>-l ,W2n-1-.l | WQ=0) . 


Each  path  from  (0,0)  to  (2n-l,-l)  must  have  probability  pnqn“*  and  by  the  ballot 
theorem  (Feller  1968,  p.66)  there  are  (2n-l)"*C(2n-l,n-l)  such  paths,  completing 
the  proof  of  (i) , 

On  the  othexhand,  each  path  from  (0,0)  to  (n,m)  has  (n*m)/2  "up**  steps 
and  (n-m)/2  "down"  steps,  thus  has  probability  p(n"ri)/2q(n+ra)/2^  gy  tjje  ballot 
theorem  the  number  of  paths  from  (0,0)  to  (n,m)  lying  above  or  on  zero  is 
((“♦l)/(n+l) }  c(n+l,(n+m+2)/2).  Thus  (ii)  follows.  || 
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